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Abstract 



We generalize the multilevel algorithm of Liischer and Weisz to 
study SU(.ZV) Yang-Mills theories with the tree-level improved gauge 
action. We test this algorithm, comparing its results with those ob- 
tained using the Wilson action, in SU(3) and SU(4) Yang-Mills theories 
in 2 + 1 and 3 + 1 dimensions. We measure the static quark potential 
and extract the Liischer term, predicted by the bosonic string theory. 

1 Introduction 

An efficient lattice gauge theory algorithm proposed by Liischer and Weisz, 
the multilevel algorithm, has been shown to provide an exponential reduc- 
tion of the statistical errors in calculations of the Polyakov loop correlation 
function p]. Multilevel algorithms are useful in many contexts, in which one 
has to cope with an exponentially decaying signal-to-noise ratio, for example 
in the computation of the glueball spectrum [2] and in the computation of 
the correlation functions, which are related to the transport coefficients of 
the quark-gluon plasma [31 S]- In [5], a multilevel algorithm was used to 
study the localization properties of gauge fields on domain wall defects. In 
this paper, we generalize the multilevel algorithm to study the static quark 
potential. 

In a confining theory, the static quark potential has the large distance 
asymptotic expansion [6] 



where a is the string tension, [i a constant (a regularization-dependent 
mass), and 7 is the Liischer term 
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with D as the dimension of the space-time [7J E] • The interaction between 
the static quark-antiquark pair can be described by an effective string theory. 
It has been suggested that the expectation values of large Wilson loops have 
a correspondence with amplitudes of an effective bosonic string theory [61 [9] . 
The Polyakov loop correlator can be used similarly to the Wilson loop to 
study string effects [Iffi ITU H2"] . 

One of the consequences of the effective string description at zero temper- 
ature is the presence of the term proportional to 1/r, i.e. the Liischer term 
Q, in the long distance inter-quark potential. The Liischer term includes a 
coefficient that depends only on the dimension of the space-time and is not 
influenced by higher order corrections of the effective string action. The aim 
of this paper is to study the potential, namely compute the Liischer term, 
in pure glue SU(iV) Yang-Mills theory, using numerical simulations. 

Although the validity of the asymptotic expansion ([!]) can be checked by 
means of numerical simulations, the problem in these lattice computations, 
is that the signal-to-noise ratio decays quickly at large distances. This makes 
it difficult to clearly separate the j/r correction from the other terms. Still, 
lattice gauge theory can offer support for the string model. For example, 
the data for the potential in a range of distances from 0.4 to 0.8 fm has been 
noted to agree with within small errors |13j . Also, using highly efficient 
simulation techniques it was confirmed [14] that the expectation values of 
large rectangular Wilson loops are matched by the string theory amplitudes 
to a precision where the subleading string effects can be observed. However, 
a more recent study [15J on string effects at large-A using Wilson loops 
seemed to differ from the effective string prediction, contrary to previous 
studies pSlin|. 

In this paper, all simulations are performed utilizing the multilevel al- 
gorithm, and besides the standard Wilson action, we will also conduct sim- 
ulations with the tree level improved action |18[ [191 [20"1 [21] . The improved 
action can strongly reduce the lattice artifacts, and hence give better re- 
sults also on relatively small lattices. Examples of "success stories" of the 
improved action for the gauge action include computations of the SU(3) 
equation of state [221 E3] and the studies of renormalized Polyakov loops 
|24[ [2"o] . Similarly, improved actions for the quarks are particularly useful 
when it comes to computationally challenging problems, like, e.g. lattice 
studies of walking technicolor models (see, for example, [20 j. 

The string-like features of the interquark 'flux tube' have already been 
studied with calculations in Z 2 [H [27J EH] , Z 4 [29j , U(l) [301 ED [32] , SU(2) 
031 EH ESI E5], SU(3) P ESI E], as well as with SU(iV > 3) [HI HJ EEl E2] , 
and even in a random percolation model (with an appropriate definition of 
the observables) |40[ 141] , By extending the research also to bigger gauge 
groups we pursue two main subjects of interest. First of all, since large- A 
studies are largely motivated by the fact that they can provide a mathe- 
matically simplified frame for studying QCD, we want to see whether the 
numerical results for the flux tube obtained with a large-A theory agree 
with SU(3) results. A second motivation comes from the conjectured con- 
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nection between large- N conformal gauge theories and string theory (see for 
example the review [12] and [13]), affiliating with the Maldacena conjecture 
and AdS/CFT. In this paper, the groups SU(3) and SU(4) are studied. 

This paper is organized as follows: In section 2, as an extended introduc- 
tion on the topic of this paper, we will discuss some basic issues on bosonic 
string theory, namely the Nambu-Goto string and how it can model the 
behavior of the flux tube in a confining gauge theory, and also some ear- 
lier lattice results relevant to this topic. Then, in section 3 we discuss the 
Polyakov loop correlation function on the lattice and the actions used in this 
work, and in section 4 we describe the multilevel algorithm in more detail. 
In sections 5 and 6 we present the results of this work and the conclusions. 



2 Bosonic string theory 
2.1 Nambu-Goto string 

The simplest effective action for a bosonic string is simply the string tension 
a times the area of the string world sheet, i.e. the Nambu-Goto action 



S e s = cr J d 2 ^detg a p. (3) 

In the so called "physical gauge" the integrand reads 

y/l + (d h) 2 + (c^) 2 + (d h x dxhf (4) 

where h is the displacement of the world sheet surface in the transverse di- 
rections. Expanding this in a perturbative series in l/(crr 2 ), at leading order 
we get the following expression for the Polyakov loop correlation function 
(see, e.g. [30J for a discussion) 

p—crrL— fiL 

( P *(r)P(0)) = - - lXxD _ 2 , (5) 



'l»2r 

where we have used Dedekind's ij function 



v(r) = q^U(l-q n ), q = e 2 ™ . (6) 



n=l 



When £- 3> 1, eq. \b\ gives the Liischer term in the quark-antiquark poten- 



2r 

tial. 



The spectrum of the Nambu-Goto string can be obtained by canonical 
quantization |47l 14*5] : the energy levels for a string with fixed ends are 



2tt / D - 2 

— k n 

ar 2 V 24 

(D-2- 24n)vr | :i 



E n (r) = ar X /l + -2(n-^^), n G N, (7) 



ar-"- — >- + (l/r% (8) 
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As a consequence, the partition function describing the string with fixed 
ends reads 

oo 

Z = J2"ne- E ^ L (9) 

n=0 

where oj n are the number of states. 

As mentioned, one of the outcomes of the effective string description at 
zero temperature is the presence of the Liischer term in the long distance 
inter-quark potential. One should not confuse this with a Coulomb term 
originating from the one-gluon exchange process. The Liischer term is a 
Casimir effect, which is due to the finiteness of the interquark distance r. As 
a heuristic argument for the term, let us consider the string-like tube created 
by two color sources separated by distance r. Due to quantum fluctuations, 
the string can vibrate, and we can express its wavelength through 

r = A(— ) (10) 

where n = 0, 1, 2, 3, . . . . If we consider each of these modes as quantum 
mechanical harmonic oscillators, then their energies are 

E k n = ho n (k+ 1 ^ (11) 
and if we just consider the ground state, k = 0, we can write 



E n 



(12) 



oo 1 1 oo oo 1 

n=0 n=0 n=0 

it fl 3 5 

= -(- + H h2H h. 

rV2 2 2 

= ^-(l + 2 + 3 + 4 + 5 + .. 

1 vr 
~2 127 

where we have used the Riemann zeta regularization. This is just for one 
transverse direction; in D dimensions there are D — 2 transverse directions, 
so in four dimensions, we have 

E n = - — . (13) 
12r V 1 

Finally, the effective string model also gives a prediction for the form of the 
inter-quark potential in finite temperature gauge theories [351 133 HH] • 

2.2 QCD string 

In a pure gauge theory the ground state interquark potential V(r) of a heavy 
QQ pair can be expressed in terms of the two-point correlation function G(r) 
of Polyakov lines 

V(r) = -ilogG(r) = ~ log(P*(r)P(0)), (14) 
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where r is the interquark distance and L the system size in the time-like 
direction [30J. 

In the confining regime of SU(iV) gauge theories, the asymptotic behavior 
of V(r) at large distances is a linear rise, and the flux lines between the 
well separated color sources are expected to be squeezed in a thin, string- 
like tube [7J El [50]. This "confining string" can be considered as a basic 
object, in particular when it cannot break, i.e. when there is no matter in 
representations of non-zero iV-ality. By studying the low-energy effective 
action on the string we pursue to understand the low-energy dynamics of 
long strings. In large- iV gauge theories the confining string can be thought 
of as a weakly coupled fundamental string subject to some effective action 
which, if known, can be used to study the low-energy properties of the 
model. Unfortunately, in general the effective string action is not known, 
but Monte Carlo computations on a lattice may offer a way to study the low 
energy effective action of a confining string and give insights of its properties. 
Furthermore, the general properties of the effective string can be derived 
based on the symmetries it should have [5H E21 E31 EH E2] • 

Even though the interquark potential V(r) and related quantities have 
been studied extensively on the lattice for many years (see, for example, 
[57} l56l [57] for references) , the question whether the picture given by effec- 
tive string models is quantitatively satisfactory, is still under investigation. 
Previous results |371I56| have shown prominent support for the bosonic string 
prediction, in particular, the Luscher term has been shown to be a universal 
feature of the IR regime of confined gauge theories. However, results from 
recent, high precision Monte Carlo simulations [5?ll5gll6ll271l531IMI50"ll53l 
l6"T] l62j suggest that higher order corrections to V(r) might not be universal. 
In fact, there are theoretical arguments suggesting that the effective string 
action is different from the Nambu-Goto one at higher orders. 

Another aspect of interest, is the excitation spectrum description. Ac- 
cording to the bosonic string picture, the excitations are expected to be 
described by the spectrum of harmonic oscillator energies 



In the ground state potential the Luscher term is supported by numerical 
evidence down to very short distances, but for the excited states the lattice 
results are not in such a good agreement with the theoretical expectations, 
as discussed in [30]. Mismatches between effective string predictions [10] 
and numerical results have been found in [631 ESI EU EE] . 

In this paper however, we shall focus the attention onto the ground state 
interquark potential. Assuming that the low energy dynamics of the pure 
gauge model is described by the effective string, and assuming that decays 
of excited states do not play a significant role, we can write the Polyakov 
loop two-point correlation function as a string partition function 



E = E + -n 



n G N. 



(15) 



r 




(16) 
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where S e g is the effective action for the world sheet spanned by the string. 

The dynamics of the confining string is not known, but it should respect 
the expected rotational symmetries. This implies that only the terms that 
are rotationally symmetric can be part of the effective string action. In 
fact these constraints are more general and they restrict the form of the 
effective action (at least at the lowest orders) to be the Nambu-Goto one. 
Lattice simulations for pure Yang-Mills theories in D = 3 and D = 4 show 
the effective action being very well approximated by this form. However, 
deviations from Nambu-Goto can be derived at higher orders (see |65j and 
references therein). 

Different approaches allow one to constrain the effective action of a con- 
fining string. Polchinski and Strominger [66, Wl\ EHj proposed to consider 
the degrees of freedom in the effective action as the embedding coordinates 
of the string in a conformal gauge world sheet. Requiring the effective action 
to have the correct critical central charge, one is left with constraints that 
have been shown to imply the four-derivative effective action to agree with 
the Nambu-Goto form. Unfortunately, generalizing this approach to higher 
orders appears to be challenging. 

Another approach [7J El |69j is to write the effective action in static 
gauge, such that the degrees of freedom are only the (D — 2) transverse 
fluctuations of the string world sheet. Essentially, the action should non- 
linearly realize the Lorentz symmetry rotating the direction in which the 
string propagates, and the transverse directions. Following the suggestion 
of Luscher and Weisz in [69| , the effective action is constrained by computing 
the partition function of long closed strings, winding around a periodic size 
of the system, and writing it in terms of a sum over string states. In [69J it 
was shown that in D = 3 the partition function on the annulus constrains 
the four-derivative terms to be of Nambu-Goto form. 

2.3 Some earlier lattice results 

In recent years the effective action of confining strings has been studied on 
the lattice with higher and higher precision (see [57] for a review of results). 
Studies of the three dimensional pure Yang-Mills theory have produced very 
accurate results of the spectrum of confining flux tubes in large- N SU(A^) 
gauge theories [381 ES]- Considering torelons of length L, it has been found 
that the ground state energy agrees with the Nambu-Goto result at order 
1/L, and is consistent with it at order 1/L 3 . Also, a possible deviation at 
order 1/L 5 comes with a very small coefficient. The excitation spectrum has 
also been found to agree with Nambu-Goto at orders 1/L and 1/L 3 , however 
at higher orders there are deviations. The lattice data are now reaching high 
precision, so that it is becoming possible to determine at which order the 
deviation occurs, which could be already at order 1/L 5 \70\ 15 lj . 

Similarly, in simulations of large- N gauge theories done in 3 + 1 dimen- 
sions, there is agreement with Nambu-Goto for large L, but the order in 
which the deviations arise is not clear. Such is the case also in simulations 
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of interfaces in the 2 + 1-dimensional Ising model, in which one can also 
see a good agreement with Nambu-Goto |71j . but still more precision would 
be needed in order to tell at what order deviations from Nambu-Goto arise 
[70]. 

2 + 1 dimensional confining strings in higher representations, called u k- 
strings", have also been studied [721 [73] & t large N. Comparison to Nambu- 
Goto showed large deviations for all states, including the ground state, pos- 
sibly starting already at order 1/L 5 . However, there are some technical 
aspects to be taken into account. In the large- N limit, assuming k to be 
fixed, the binding energy of k fundamental strings to form a A;-string may 
vanish as 1/N 2 or as 1/N |74j . implying that in the large- iV limit there are 
at least (k — 1)(D — 2) light modes on the worldvolume of a fc-string, whose 
mass goes to zero in the large- N limit as 1/N or as 1/N 1 / 2 , respectively. 
The relevant length scales concerning the form of the effective action are 
larger than N/\/T (or N 1 / 2 /\/T). 

Finally, we mention that the ground state potential of a confining string 
in the continuum limit of a percolation model in 2+1 dimensions was studied 
in [IT], and again agreement was found with Nambu-Goto at orders 1/L and 
1/L 3 , however deviations were seen at order 1/L 5 . As discussed in |70] , this 
model is interesting in itself, as it is not necessarily expected to correspond 
to a weakly coupled string theory. 



3 Polyakov loop correlation function on the lattice 

We discretize the SU(iV) Yang-Mills theory on an isotropic 4-dimensional, 
or alternatively, 3-dimensional lattice with spacing a, time-like extent aNt 
and spatial size aN s . Apart from providing precise results more quickly, 
simulations in 3 dimensions are interesting because there the 1/r term is 
certain to be the Liischer term; in 4L> it could get mixed with a Coulomb 
interaction. In 3D the Coulomb interaction is of the form log(r), rather 
than 1/r. 

In all directions of the lattice we impose periodic boundary conditions. 
For the gauge field action we take the usual Wilson plaquette action 

S = ' 3 E (l-^ReTrLV(x)), (17) 

where U^ u is the plaquette oriented along the ^-plane and located at the 
lattice site x, and j3 = 2N/g 2 . We also carry out the simulations using the 
tree-level improved action [TBI [T9| 1201 [21] 



S = P E I 1 " ^ ReTr \l U l^ x ) - ^ 2 (x) - ±U$(z)\ | 

.I'.//' v L J ' 



(18) 



For any given gauge field configuration U (x, fj,), the (trace of the) Polyakov 
loop is defined as 

P(x) = TV U(x, n)XJ{x + an, n) . . . U(x + (T - o)£, /x) M=0 (19) 
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and the correlation function 

<P(*)*P(l/)> = | f HdU(x^)P(xyP(y)e- s M. (20) 

The ground state inter-quark potential V(r) of a heavy QQ pair in a pure 
gauge theory can be expressed with the two-point correlators of the Polyakov 
loop [H] 

V(r) = -- i -log(P(xyP(y)) + e, (21) 
aN t 

where 

e = ^-(oj ie - AET + ...), AE = E 1 -E . (22) 



4 The Multilevel algorithm 

Like mentioned in the introduction, the lattice data is obtained using the 
multilevel algorithm of Liischer and Weisz p]. In this algorithm, using 
the locality of the theory, the lattice is split into sublattices that do not 
communicate with one another, and the final observables are built combining 
the independent measurements of each sublattice. Of course, from time to 
time the boundaries between the sublattices are updated, so that the final 
results are the same as in the usual theory. 

Following the discussion in [33] , suppose we measure an observable O by 
combining the results of averages O su b computed in N different sublattices. 
With J\f sublattice measurements, the combination of the sublattice averages 
O su b corresponds to (N)^ measurements of O, i.e. we get an estimate of 
O as if (N)-^ measurements would have been performed. However, due to 
the links that have been kept frozen at the boundaries of the sublattices, 
this estimate is biased by a background field, but this bias is removed by 
averaging over the boundary configurations. 

Using the notation of [33J, consider the correlation function of two Polyakov 
loops 



^(0,0),4 • • • "(5,1,-1) ,4 



(P(0)P(x)*> = i [ HdUy^Tr 

[uko)A---UkL-i) A y S[U] - (23) 



x TV 



Next, we slice the lattice along the temporal direction into N = Nt/nt 
sublattices, i.e. nt is the temporal thickness of each sublattice in units of the 
lattice spacing a. Now, to obtain sublattices isolated from each other, we fix 
the set V{! of all spatial links with time coordinates krit, k = 0, . . . , (Af — 1). 
This way, the dynamics of every sublattice depends only on the background 
field of the two frozen time slices, and hence are totally independent from 
one another. 
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As in [33] , we rewrite ( 23 ) as 



(P(O)F(x)*) 



n 



dui s) T-, 



6,(x) [V o 



Mi 



where 



rWl 



(24) 



rpasyPSryis) y(s)i 

8,(x) [ « ' J J 



n ^ 



^(0,0) ,4 ' ' • ^(0,n t -l),4 



^(*x,0),4 ' 



07 



17, 



(x,n t -l),4 



e -S[U;V^,V^] 

V Z[V} 8) ,V} 3) ] ' 



(25) 



The partition function of the sublattice with fixed temporal boundaries 
reads 



zwr,v 



in which S[U; V i , Vj ] is the action in the sublattice with fixed temporal 



(26) 



y,v 



boundaries and V} S K a, f3, 7 and S are color indices and T" 7 ' 35 are 



gauge-invariant quantities under sublattice gauge transformations. V[V k 
is the probability for the spatial links with time coordinates knt, k 
0, . . . , (Af — 1), to be Vl , and can be written 

nvt s) ] =u n^n^ - v k (s) )e- s[u] 

J y,A» k 



Ml 



(27) 



Iterating this procedure and averaging the background field according to 



the probability distribution of eq. (27) gives a numerical estimate for the 



integral (24) [33J. This technique is the so-called single level algorithm. A 



more general "multilevel algorithm by Liischer and Weisz involves a feature, 

(s) 

that the updating frequency of the background field V fc v is not the same 
for the various time slices k. However, for many purposes the single level 
algorithm seems to be more efficient E3] • 

In numerical simulations with the single level algorithm, three parame- 
ters have to be fixed: the temporal sublattice thickness nt, the number TV" 
of sublattice measurements and the number M of background field config- 
urations to integrate over V^ s \ These parameters are dependent on each 
other as well as on the temporal extension Nt of the whole lattice and on 
the distance between the two Polyakov loops. Finding the optimal choice 
for them is thus not trivial. Some results on the optimization step can be 
found in [MIE] • 

When measuring the Polyakov loop correlation function (P(O)P(x)*) 
the final error bar is the combination of the uncertainties of the sublattice 
averages and their fluctuations, due to different background fields. With nt 
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fixed, a large distance between two loops requires both N and M to be large. 
N is typically order of several thousands and M of few hundreds. It should 
be noted that N does not depend on Nt whereas M does. With the Liischer 
and Weisz algorithm, exponential gain in the accuracy of the numerical 
estimation of (P(O)P(x)*) is possible only in the temporal direction. Every 
sublattice estimate decreases exponentially with the distance, but it is still 
estimated with an error reduction proportional to 1/ t/N [33J . 



Wilson action 



Improved action 



N t 



N, 



/ 7 



Figure 1: For an example on how the multilevel algorithm works, i.e. which 
links are updated during the sublattice updates, consider a lattice divided 
into three sublattices. On the left, a presentation of the Wilson action, 
and on the right, the improved action. On the borders of a sublattice, the 
red links are the ones kept fixed, whereas the blue links are updated. The 
improved action is slightly less local than the Wilson action, because it 
involves not only plaquettes but also rectangles. However, like the Wilson 
action, the improved action is still "ultralocal" , namely, in the action, each 
link is coupled directly to other links only up to a finite distance. 



5 Results 

The simulations were performed with gauge groups SU(3) and SU(4), in 



2 + 1 and 3 + 1 dimensions, with the general Wilson action (17) and with 



the tree- level improved action (18). The totalling 8 different cases (up to 4 
different values of (3 for each) were studied using the multilevel algorithm 
with sublattice size of 4 lattice spacings. In 3D with SU(3) the lattice size 
was N t = 48, N s = 24, in the other cases N t = 24, N s = 16. The first 1000 
updates were used for thermalization, and the following 1000 — 5000 updates 
were used for the actual measurements of the Polyakov correlation function. 

In addition to considering the correlators measured with an "on-axis" 
increasing distance between the sources, we also measure the correlators in 
"off-axis" points, i.e. where the sources have "vertical" distance in addition. 
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In general, the lattice breaks rotational symmetry, so off-axis correlators 
of Polyakov loops are not expected to lie exactly aligned with the on-axis 
ones. The full continuum rotational symmetry is however restored when the 
lattice spacing tends to zero [75J. In the off-axis points one could expect 
to see a difference between the results obtained using the general Wilson 
action and the tree-level improved action; when calculating the potential, 
the off-axis points are expected to fall further from the expected curve with 
the non-improved action. This makes measuring the off-axis points a valid 
test to assess the "goodness" of the improved action. 

For each value of /3, the static quark potential was calculated according 
to (21), see figures 2|9 In 4D cases, a function of the form 



V(r) = ar + V + j/r (28) 

was fitted to the data, thus giving estimates for a, Vq, and 7. 

In 3D cases the fitted function could basically include also a logarithmic 
term, i.e. the Coulomb term, which would play a role on small distances. 
However, including this term makes the fit very unstable on larger distances, 
as seen also in [76]. Instead, following the conclusions of [69], in 3D cases 
we add a 1 /r 3 term in the fitting function, 

V(r) = ar + V + j/r + b/r 3 . (29) 



Compared to (28), this form gives us better a fit and better estimates for 
the Liischer term 7. 

In 3-D the bosonic string prediction for the Liischer term is 7 ~ 0.1309, 
and in AD 7 ~ 0.2618. As we can see from the tables [TJ [2j in both 3D 
and 4D, the Liischer term extracted from our lattice data follows these 
predictions. 

All statistical uncertainties affecting the results were computed using 
standard statistical analysis methods. In particular, the error bars on the 
values of the average correlators were computed with the jackknife method, 
while those on the potential were calculated by Gaussian propagation of the 
latterH 

5.1 Systematic uncertainties 

As always with results from lattice studies, there are some systematic un- 
certainties that need to be taken in account. Finite volume effects, finite 
cut-off effects and systematic effects in the choice of the fitting function are 
all sources of such uncertainties in our calculations. 



lr This procedure turns out to be reliable for our data; computing the uncertainties on 
the potential using the jackknife method leads to comparable results: the small differences 
between the results obtained with the two methods do not have a really significant impact 
on the findings of the study, and we include them among the possible sources of systematic 



uncertainties that we discuss in subsection 5.1 
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The finite volume effects are exponentially suppressed as e~ maL , where 
mc is the mass of the lightest glueball, and L the shortest lattice size. e~ m ° L 
is already very small when mcL > 4, so what we want is 

L>— • (30) 

The lattice sizes chosen in this work meet this requirement. 

As usual, the finite cut-off effects can be controlled by performing simu- 
lations at several values of the lattice spacing. The calculations in this paper 
involve a few different values of /3 in each studied case, and the results indi- 
cate the impact from cut-off effects to be very small. As can be seen from 



figure [TOj the results are clearly already quite close to the continuum, be- 
cause the points obtained from simulations at different lattice spacing fall on 
the same curve and because there is no significant difference between points 
obtained from on-axis versus off-axis correlators (hence the full rotational 
symmetry is approximately restored). 

The systematic effects in the choice of the fitting function can be studied 
by extending it with an extra term, c/r a with a = 2 or 3, and seeing how 
the fit parameters behave. With SU(3) in 3D and 4D, adding the extra 
term with a = 2 changes aa 2 and V$a only slightly (for example in 3-D 
with (3 = 18.0, aa 2 changes from 0.04482(40) to 0.0470(59) and V a from 
0.2255(23) to 0.202(61)), but gives very poor estimates for 7. With SU(3) 
in AD adding the extra term with a = 3 gives poor estimates for all the 
parameters; with j3 = 5.9, aa 2 changes from 0.0904(70) to 0.08(13), V$a 
from 0.670(21) to 0.72(69), and 7 from -0.236(14) to -0.31(97). 

The quality of the fit and, more precisely, the x 2 are heavily influenced by 
the choice of points that are included in the fitting process. Small distance 
off-axis points have a huge impact; for example with SU(3), (3 = 18.5 in 3D, 
we get x 2 = 1-04 when the small distance off- axis points are neglected in the 
fitting, but when included, the x 2 grows to a staggering 4445.21. Obviously, 
this is due to systematic effects related to lattice artifacts at short distances. 
Already removing just a couple of these points brings the x 2 down to 107.028. 
In another example, with SU(4), j3 = 23.0 in 3-D, we get x 2 = 0.91 with 
the small distance off-axis points neglected, and x 2 = 3627.49 when they 
are included. For reliable results, all fits in this paper are done excluding 
the small distance off-axis points. The points are however included in the 
figures 



6 Conclusions 

In this work, we have generalized the Luscher and Weisz multilevel algorithm 
for the tree-level improved gauge action to study SU(-/V) Yang-Mills theo- 
ries. Applying this efficient lattice gauge theory algorithm, we computed 
the static quark potential in SU(3) and SU(4) Yang-Mills theories in 2 + 1 
and 3 + 1 dimensions. The calculated Luscher term was shown to satisfy 
the predictions of bosonic string theory, more or less equally well with the 
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D 





aa 2 


V a 


7 
/ 


b 


3 


3 


18.0 
18.5 
19.0 


0.043927(57) 
0.041339(58) 
0.039042(62) 


0.23041(32) 
0.22689(33) 
0.22313(36) 


-0.11434(47) 
-0.11276(48) 
-0.11069(53) 


0.02313(23) 
0.02306(21) 
0.02276(23) 




4 


5.9 
6.0 


0.0904(70) 
0.0649(21) 


0.670(21) 
0.6788(63) 


-0.236(14) 
-0.2412(42) 




4 


3 


32.0 
33.5 
34.0 
30.0 


0.04716(15) 
0.04244(25) 
0.04106(27) 
0.054518(76) 


0.31280(90) 
0.3070(15) 
0.3050(16) 
0.32240(44) 


-0.1230(13) 
-0.1210(23) 
-0.1202(25) 
-0.12795(64) 


0.02528(60) 
0.0255(10) 
0.0254(11) 
0.02581(28) 




4 


11.0 


0.0680(14) 


0.737(19) 


-0.255(12) 





Table 1: fit results with Wilson action 
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aa 2 
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b 


3 


3 


18.0 
18.5 
19.0 


0.034971(85) 

0.03320(14) 

0.031427(70) 


0.20028(48) 
0.19695(76) 
0.19493(40) 


-0.10488(68) 

-0.1025(11) 

-0.10222(58) 


0.02115(29) 
0.02061(48) 
0.02106(25) 




4 


4.5 


0.0536(22) 


0.6686(66) 


-0.2478(44) 




4 


3 


23.0 
23.5 
24.0 


0.07118(38) 
0.06852(47) 
0.06599(24) 


0.3191(21) 
0.3145(26) 
0.3102(13) 


-0.1333(31) 
-0.1292(38) 
-0.1250(19) 


0.0244(13) 
0.0232(17) 
0.02187(85) 




4 


8.0 


0.058(13) 


0.771(41) 


-0.290(28) 





Table 2: fit results with the improved action 



Wilson action Improved action 



N c 
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P 


reduced \ 2 


P 


reduced \ 2 


3 


3 


18.0 


1.53 


18.0 


0.26 






18.5 


2.25 


18.5 


1.04 






19.0 


3.48 


19.0 


0.42 




4 


5.9 
6.0 


0.53 
1.21 


4.5 


0.49 


4 


3 


32.0 


15.91 


23.0 


0.91 






33.5 


65.21 


23.5 


2.07 






34.0 


86.83 


24.0 


0.79 






30.0 


1.81 






4 


11.0 


0.55 


8.0 


0.98 



Table 3: The reduced x 2 f° r each of the fits. 
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Figure 3: SU(3), in 3 dimensions, with the improved action. 4000 measure- 
ments, lattice size 48 x 24 2 . 
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Figure 5: SU(3), in 4 dimensions, with the improved action. 3000 measure- 
ments, lattice size 24 x 16 3 . 
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Figure 7: SU(4), in 3 dimensions, with the improved action. 5000 measure- 
ments, lattice size 24 x 16 2 . 
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Figure 9: SU(4), in 4 dimensions, with the improved action. 1100 measure- 
ments, lattice size 24 x 16 3 . 
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Figure 10: The physical values calculated from the results obtained with 
SU(3) in 3 dimensions, with the improved action 



Wilson action as with the improved action. Also, we can see that the results 
obtained with SU(4) agree with the results of SU(3), and in fact, the bosonic 
string prediction for the Liischer term is met more accurately with SU(4). 

Besides the calculations of the static quark potential, multilevel algo- 
rithms can be used in many other contexts, such as in the computation of 
glueball masses and the correlation functions related to the transport coef- 
ficients of the quark-gluon plasma. For future endeavors, in principle the 
multilevel algorithm for the improved gauge action could also be used in 
graphics processing unit (GPU) implementations. As an example of this 
type of application, very recently this has been done for the compact U(l) 
lattice gauge theory with the standard Wilson action, in [77] , 

Another possible useful application for a high-precision multilevel algo- 
rithm with the improvement in 3-D could be in the context of dimensionally- 
reduced effective theories for hot QCD, which, for example, very recently 
have been used to compute certain non-perturbative contributions to the 
jet-quenching parameter in [75} 175]. 
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